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ABSTRACT 

We present time-resolved spectroscopy of the eclipsing, short period cataclysmic 
variable CTCV J1300-3052. Using absorption features from the secondary star, we 
determine the radial velocity semi-amplitude of the secondary star to be K% = 378 
± 6 km s" 1 , and its projected rotational velocity to be wsini = 125 ± 7 km s _1 . 
Using these parameters and Monte Carlo techniques, we obtain masses of M\ = 0.79 
± 0.05 Mq for the white dwarf primary and M 2 = 0.198 ± 0.029 M for the M- 
type secondary star. These parameters are found to be in excellent agreement with 
previous mass determinations found via photometric fitting techniques, supporting the 
accuracy and validity of photometric mass determinations in short period CVs. 

Key words: binaries: close - binaries: eclipsing - stars: dwarf novae - stars: low mass, 
stars: novae, cataclysmic variables - stars: evolution. 
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1 INTRODUCTION 

Cataclysmic variable stars (CVs) are a class of interacting 
binary system undergoing mass transfer from a Roche-lobe 
filling secondary to a white dwarf primary, usually via a gas 
stream and accretion disc. Through determinations of the 
masses and radii of the component stars in CVs, it is possible 
to test fundamental th eories regarding their formation, ori- 
gin a nd evolution (e.g. iLittlefair et al]|2008l ; ISavourv et al.l 

OM: 

I Savoury et al.l (|201 lh carried out a photometric study of 
eclipsing CVs and found the masses and radii for both the 
white dwarf and donor star in 14 systems. These masses were 
found by fitting a parameterised model to the eclipse light 
cur ves. Th i s mod e l is ba s ed on t he techniques develope d 
by iBailevI l l 19791) ISmakl dl979h. ICook fc Warner! JjWM ), 



Wood et all (ll985l1.IWood et all lfl98^1.lHorne et all Jl994h . 
Littlefair et al.l l|200ct ) and Copperwheat et al.l ( 201(1 ) and 



relies on just four assumptions: the bright-spot lies on the 
ballistic trajectory from the donor star; the donor star fills 
its Roche lobe; the white dwarf is accurately described by 
a theoretical mass-radius relation; and the whole of the 
white dwarf is visible with an unmodified surface bright- 
ness. It is photometric mass and radii determinations such as 
these that are used to calibrate the t-q ( super-hump ex cess- 
mass ratio) relation s of lPattersonl (|2005t ).lKnigge (200fJ) and 
iKnigge et al](|201ll ), which can then be used to derive donor 
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mass estimates for large samples of CVs. It is therefore im- 
portant to check the validity of photometric mass determi- 
nations. 

For objects with periods above the period gap (the 
dearth of systems between 2.2 and 3.2 hours, see e.g. 
Kolb & Ritter 2003, Knigge 2006), the photometric fit- 
ting technique appears robust, with donor star radial ve- 
locities predicted by phot ometry in agreement with those 
found by other technique s (|Watson et al.ll2003l ; lFelinell2005l ; 
ICopperwheat et al.ll2010T ). However, for objects below the 
period g ap, independe n t test s of the photometric technique 
are rare lTulloch et all j2009) found the radial velocity of 
the white dwarf (Ki) in SDSS J143317. 78+101123.3 (P orb 
= 78.1 mins), as measured from disc emission lines, to be in 
excellent agreement w ith the photometric value predicted by 
ILittlefair et all (|2008l U The agreement is encouraging, but 
the motion of the inner disc does not necessarily follow the 
motion of the white dwarf, and so K\ estim ates from dis c 
emission shoul d be treated with caution (e.g. lMarsbJ [l988'). 
More recently, ICopperwheat et al.l |201ll 'l found the radial 
velocity and rotational broadening of the secondary star in 
OY Car (P or b = 90.9 mins) to be in good agreement with 
those predicted by photo metric methods |Wood fc Hornd 
1 199d : ILittlefair et al. 2008). However, such is the importance 
of mass determinations in CVs, additional verification across 
a range of orbital periods is highly desirable. 

One of the systems observed by I Savoury" et all (|201lh 
was CTCV J1300-3052 (hereafter CTCV 1300). CTCV 1300 
is a dwarf nova that was discovered as part of the Calan- 
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Tololo Survey follow up l|Tappert et al.ll2004I ). It was found 
to be eclipsing, with an orbital period of 128.1 minutes, plac- 
ing it immediately below the period gap. The average spec- 
trum showed clear emission lines from the accretion disc and 
absorption lines from the donor star. It is through absorption 
lines such as these that we can determine the radial velocity 
and rotational broadening of the secondary star, which can 
in turn be used to derive an independent me asure of the 
masses and radii of the component stars (e.g. | Horne et al.l 
1 19931 ; ISmith et al.lll998l : iThoroughgood et al.ll200ll . I2004T ). 

In this paper we present time resolved-spectroscopy of 
CTCV 1300 and determine the system parameters. The pa- 
rameters derived using spectroscopy will pro vide an indepen- 
dent t est of the photometr i c meth ods used bv lLittlefair et all 
(|2008h and lSavourv et al l (|201lh . 



2 OBSERVATIONS 



CTC V 1300 was observed using X-shooter (|D'Odorico et all 
2006) in service mode mounted on UT2 (Kueyen) on the 
8.2-m Very Large Telescope (VLT) on the nights beginning 
9 Feb 2010 and 6 March 2010. In total, we obtained 48 spec- 
tra (24 on each night) covering 1.5 orbital cycles, and a 
wavelength range of ~3000-24800 A. Exposure times were 
235 seconds in the UVB-arm (3000-5500 A), 210 seconds in 
the VIS-arm (5500-10000 A), and 255 seconds in the NIR- 
arm (10000-24800 A), with dead times between exposures of 
approximately 8, 9 and 1 seconds, respectively. 

The target was observed with the 1.0" xll" slit in 
the UVB-arm, the 1.2"xll" slit in the VIS-arm, and the 
0.9"xll" slit in the NIR-arm. The resolving power was 
~5100 (59 km s _1 ) in the UVB and NIR-arms, and ~6700 
(45 km s _1 ) in the VIS-arm. Seeing conditions on both 
nights were fair, varying between 0.5 and 1.5 arcseconds, 
but with flares of up to 2.0 arcseconds. 

Observations of the standard star CD 153 were used to 
flux calibrate the data and correct for telluric absorption. 
The data were obtained in 'stare' mode rather than nod- 
ding along the slit as is normal for long slit infra-red spec- 
troscopy. Consequently, the sky subtraction on the NIR-arm 
spectra is significantly worse than usual with X-shooter. We 
also obtained spectra of the spectral type templates GJ2066 
(M2V) and GJ1156 (M5V), although these data were taken 
on the nights of 11 Dec 2009 and 28 Jan 2010, respectively. 



3 DATA REDUCTION & ANALYSIS 

Data reduction was carried out using the X-shooter pipeline 
(version 1.2.2) recipes within ESORex, the ESO reduction 
execution tool. The data for all three arms were reduced 
with similar procedures. The required calibration frames 
were constructed using the standard recipes provided in the 
pipeline. In brief, they include a map of bad pixels, a mas- 
ter bias (for the UVB and VIS-arms), a master dark (for 
the NIR-arm, as dark contribution is negligible in UVB and 
VIS) and a master flat. The data was first bias and dark 
subtracted, before an inter-order background was fitted and 
subtracted. Science frames were then divided by the flat 
field, and then the object was localised on the slit. Sky sub- 
traction and cosmic ray removal took place, and the data 



then underw ent an optimal extraction routine (|Horndll986l ; 
lMarshlll989l ). and order merging. 

Wavelength calibration was undertaken using an arc 
line spectrum which was taken from the ESO archive. The 
wavelength calibration recipe used a physical model to gen- 
erate a best-guess solution of the line positions on the cali- 
bration frame. These lines were then fitted by 2D Gaussians, 
and the resulting positions of these lines were adjusted via a 
polynomial fit to the whole CCD. From residuals to the line 
fitting, we estimate that this calibration is accurate to ~1 
km s _1 (at A = 8183 A). We corrected for flexure in the VIS 
arm by measuring the shift of t he observed sky li nes relative 
to their positions measured by iHanuschikl {2003). The indi- 
vidual spectra were then moved by these shifts, which were 
typically between 10 and 35 km s~ . We do not attempt to 
correct for flexure in the UVB and NIR arms, since we do 
not attempt to measure radial velocities to a high degree of 
precision in these bands. 

The time and wavelength axis of the data were corrected 
to the heliocentre. 



4 RESULTS 

4.1 Average Spectra 

The average spectra of CTCV 1300 are shown in Fig.Q] The 
upper panel shows the wavelength range 3200-5500 A (from 
the UVB-arm), the centre panel shows 5750-10000 A (from 
the VIS-arm), and the lower panel 10000-13500 A (from the 
NIR-arm). Each spectrum is flux calibrated, and telluric cor- 
rection has been attempted. 

Throughout the spectrum we see strong, broad, double- 
peaked Balmer lines and several double-peaked He I lines 
(4471, 4922, 5015, 5875, 6678, 7065 and 10830 A). Broad 
double-peaked lines such as th ese are typical of a high- 
inclination accreting binary (e.g. iHorne fe MarsrJll986h . 

The high ionisation line He II 4686 A appears absent 
in the average spectrum, but is visible in the trailed spec- 
tra (see Sections 14.21 and I4.3f) . Several absorption lines are 
present between 4000-4800 A (see Fig. [2} which appear to 
trace the motion of the disc (see Section |4. 2 1 . We believe the 
most likely cause of these absorption lines is a veil of disc 
material along the line of sight, The majority of these lines 
appear to be Fe I, Fe II and Ca I. Simila r features have been 
obs erved in the spectrum of O Y Car by IHorne" et all (|l994h 
and ICopperwheat et al.l (|201lh . 

The helium lines at 4922, 5015 A appear to show strong, 
narrow absorption cores that dip below the continuum, as do 
the higher-order Balmer lines between 3600-4000 A. The O I 
triplet at 7773 A is clearly visible, and also appears to drop 
below the continuum. F eatures such as these are observed 
in a number of CVs (e.g |Marsh|[l987l ; IWade fe Hornel fl988; 
iFriend et al.lll988l ). These absorption cores are believed to 
originate through self-absorption in the accretion disc. 

The Ca II triplet at 8498, 8542 and 8662 A (here- 
after 8567 A) is clearly present and originates from the 
disc, although there is evidence of emission from the irra- 
diated side of the donor (see Sections 14.21 and 14. 3[) . Sim- 
ilar features have been observe d in the spectrum of GW 
Lib (Ivan Spaandonk et al.ll2~010l ). The higher orders of the 
Paschen series are also visible from ~8800 A onwards, and 
are possibly blended with the Ca II emission. 
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Figure 1. The average spectra of CTCV 1300, in the rest frame of the binary. The upper panel shows the UVB-arm, the centre panel 
the VIS-arm, and the lower panel the NIR-arm. The most prominent features are labelled. 



Absorption features from the secondary star are clearly 
visible in the form of TiO bands around 7100 A and 7600 A, 
and weak K I absorption doublet at 7664, 7699 (hereafter 
7682 A), 11773, 12432 and 12522 A. However, these regions 
are heavily affected by telluric absorption. The clearest fea- 
tures from the secondary star are the Na I doublets at 8183, 
8194 A (hereafter 8189 A) and 11381, 11404 A (hereafter 
11393 A), although the second of these is also heavily af- 
fected by telluric absorption. 



4.2 Trailed Spectra 

The data we re phase binned into 30 bins, according to the 
ephemeris of lSavourv et al.l (|201lf l. The UVB-arm has com- 
plete phase coverage although, due to the differing exposure 



times (see Section [2]), one of these bins is empty in both the 
VIS and NIR-arms. 

We divided the continuum by a polynomial and re- 
binned the spectra onto a constant velocity-interval scale 
centred on the rest wavelength of the lines. Fig. [3] shows the 
trailed spectra of the Ha, H/3, H7 and H<5 lines in CTCV 
1300. Each line shows two clear peaks that vary sinusoidally 
with phase, in addition to the characteristic s-wave between 
phases 0.1-0.4 from the bright spot. 

In Fig. [4] we show the trailed spectra of two Na I dou- 
blets (8189 and 11393 A), the Ca II triplet (8567 A), the 
7682 A K I doublet and He II (4686 A). The phases at which 
the Na I and K I lines show maximum red-shift (0 = 0.25) 
and blue-shift ((/> = 0.75) suggest they originate from the 
donor star. We see evidence for emission from the donor 
star in the Ca II lines through a component in the trail that 
is in phase with the Na I lines. However, this component 
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Velochy (km/s) Velocity (km/s) 

Figure 3. The trailed spectra of the Ha (top left), H/3 (top right), H7 (bottom left) and H<5 (bottom right) lines in CTCV 1300. 



is only visible during phases ~0. 25-0. 75, which indicates it 
arises from the irradiated side of the donor. The He II line 
appears to follow the motion of the bright spot, as defined 
by the s-wave in the Balmer trails. 

In Fig. [5] we show the trailed spectra of the absorption 



line forest between 4130-4300 A. The lines all appear to 
move together, suggesting a common place of origin. Using 
the same method outlined in Section l4T4l we find the velocity 
of these lines to be K a bs = 116 ± 4 km s" 1 , with a phase 
offset of A</> = 0.072 ± 0.006. The high velocity (compared 
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Figure 4. The trailed spectra of the 8189 and 11393 A Na I doublets (upper left and upper centre, respectively), the 7682 A K I doublet 
(top right), the Ca II triplet (8498, 8542 A bottom left, 8662 A bottom centre) and He II (4686 A, bottom right) in CTCV 1300. Black 
and white lines represent absorption and emission, respectively. 
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Figure 5. Trailed spectra of the forest of Fel, Fell and Cal ab- 
sorption lines between 4130-4300 A. 

to the expected motion of the white dwarf, ~ 90 km s _1 , see 
Section [4. 6p and significant phase offset suggests that these 
lines originate in the disc. 

4.3 Doppler tomography 

Doppler tomography is an indirect imaging technique which 
can be used to determine the velocity-space distribution of 
the emission in cataclysmic variables. For a comprehensive 



revi ew of Dopple r tomography, see iMarsh fc Hornel (|l988l ) 
and lMarshl l|200ll ). 

Fig. [6] shows Doppler maps for Ca II (8498, 8542 fc 
8662 A), Ha, H/3 and He II (4686 A). Eclipse data (between 
phases 0.95 and 1.05) are removed. A systemic velocity of 7 
= -20 km s _1 was applied to shift the maps onto the K x — 
km s _1 axis (see Section r4.4[) . 

In each map we see a ring-like distribution of emission 
centred on the white dwarf, which is characteristic of an ac- 
cretion disc. In the Ca II maps, we see an enhanced area 
of emission at velocities intermediate to the free-fall veloc- 
ity of the gas stream (lower stream) and the velocity of the 
disc along the gas stream (Keplerian velocity, upper stream) . 
This emission is attributed to the bright spot. The three Ca 
II maps all show clear emission from the donor star, which 
appears to be concentrated towards the inner hemisphere, 
indicating that irradiation is significant. The Ha map shows 
emission from the secondary star, a feature uncommon in 
short period CVs. The H/3 map shows weak bright spot emis- 
sion. The He II emission appears to show emission near both 
the Keplerian velocity stream and at velocities intermedi- 
ate to the Keplerian velocity stream and free-fall velocity 
stream, although it is possible that this is an artifac t arising 
from limited phase coverage (Ma rsh fc Hornelll988l ). If this 
is a genuine feature, its position relative to the Ca II emis- 
sion suggests that the He II emission is caused by a mixture 
of gas-stream and disc material. He II bright spot emission 
at Keplerian disc velocities has been observed in other sho rt 
period CVs |Marsh et al.lll990l ; ICopperwheat et alj|201ll ). 
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Figure 6. Doppler maps of CTCV 1300 in Ca II (8498, 8542, 8662 A), Ha, H/3 and He II (4686 A) computed from the trailed spectra 
in Figs. 151 andl4l Data taken during eclipse have been ommited from the fit. The predicted position of the secondary star and the path 
of the gas stream are marked. The three crosses on the map are, from top to bottom, the centre of mass of the secondary star, the 
centre ofmass of the system, and the white dwarf. These crosses, the Roche lobe of the secondary, the Keplerian velocity along the gas 
stream (top curve), and the predicted trajectory of the gas-stream (bottom curve) have been plotted using the system parameters found 
m Section H21 The series of circles along the gas stream mark the distance from the white dwarf at intervals of O.lLi, where l.OLi is 
the secondary star. 



4.4 Radial velocity of the secondary star 

The secondary star in CTCV 1300 is visible through weak 
absorption lines. The strongest of these lines is the Na I 
doublet at 8189 A. In order to determine the radial velocity, 
we cross-correlated the individual spectra of CTCV 1300 
against an average spectra of CTCV 1300 using an iterative 
technique. We chose the Na I line at 8189 A because it is 
much stronger and less affected by telluric absorption than 
the line at 11393 A. 

We subtracted fits to the continuum from the individ- 
ual spectra and then corrected for the orbital motion of the 
secondary star with a first guess of K 2. For each individ- 
ual spectrum, we then created a template spectrum that 
consisted of an av erage of all the spectra minus th e spec- 
trum under study (|Marsh. Robinson fc Woodlll994l ). These 
template spectra were then cross-correlated against the un- 
corrected data. The velocity shifts as a function of orbital 
phase were then fit with a sine function according to; 

V = 7-Jfafflnphr(0-0o)l (1) 

where V is the velocity shift, 7 is the systemic velocity of the 
system, K2 is the radial velocity of secondary star, <f> is the 
orbital phase, and <f)o is the phase offset. This then yielded 



a new value of K2 and 7 to correct our spectra with. We 
added an intrinsic error in quadrature to each error bar to 
account for systematic error, and reach a reduced- y 2 of 1. 
This process was repeated until K2 converged. We arrive at 
a value of K2 = 379 ± 6 km s _1 , with an intrinsic error of 
22 km s _1 added in quadrature to each error bar. The radial 
velocity curve obtained using this technique is shown in the 
upper panel of Fig. [7] 

The value of 7 obtained via our auto-correlation tech- 
nique is not representative of the true systemic velocity, 
which must be determined through cross-correlation with 
a template star of known radial velocity. Therefore, in order 
to verify this result and find the true systemic velocity (7), 
we then cross-correlated against our M-dwarf template spec- 
tra using the same wavelength range. The template spectra 
were artificially broadened by 46 km s^ 1 to account for the 
orbital smearing of CTCV 1300 through the 210-second VIS- 
arm exposures, and then by the best-fitting values for the 
rotational velocity of the secondary star (vsini) found in 
Section 14.51 An intrinsic error of 22 km s~ x was added to 
each error bar from the M2V cross-correlation, and 24 km 
s _1 to M5V data, to account for systematic errors and reach 
a reduced-x 2 of 1. The radial velocity curves are shown in 
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Figure 7. The radial velocity curve of CTCV 1300 obtained 
through auto-correlation (upper panel), cross-correlation against 
an M2 template (centre panel) and cross-correlation against an 
M5 template (lower panel). 

the centre panel (M2), and bottom panel (M5) of Fig. 
Cross-correlating against the M2 and M5 templates yield 
values of Ki = 373 ± 6 km s _1 and K 2 = 376 ± 7 km 
s -1 , respectively. The M2 template could not be corrected 
for flexure, so we only use the M5 template to derive the 
systemic velocit y for CTCV 1300 . Using the radial veloci- 
ties provided bv icizis et"aH (2002), we find 7 = -20 ± 5 km 
s _1 . For K2, we prefer the value found through the auto- 
correlation, that is K2 = 379 ± 6 km s _1 . This is because 
the average spectra of the data is a better match to the data 
than the M5 and M2 templates. 

The radial velocity curves produced through this tech- 
nique show some variation from a sine fit between phases 0.4 
to 0.6, which is characteristic of ir r adiation suppressed ab - 
sorption (e.g. lBillington et al.|[l99rj ). lMarsh fc Home! (|l988l ) 
recommend only fitting the above data between phases 0.8 
to 1.2, since at these phases the effects of irradiation are at 
a minimum. Fitting the auto-correlation data, we obtain a 
value of K2 — 378 ± 6 km s~ , which we use hereafter. This 
value is consistent with the value predicted by I Savoury et all 
|201ll ). K 2 = 372.2 ± 2.5 km s _1 . Fitting the M5 and M2 
templates between the same phases gives K2 = 372 ± 7 km 
s _1 and K2 = 378 ± 8 km s _1 respectively, which is again 
in excellent agreement with the auto-correlation data, and 
consistent with the photometric method. 

4.5 Rotational velocity of the secondary star 

The normalised spectra of CTCV 1300 were corrected for 
the orbital motion of the secondary star using the value of 
K2 obtained in Section T4. 41 The spectra were then averaged 
together in order to maximise the strength of the Na I dou- 



blet at 8189 A. The spectral-type templates were broadened 
to match the smearing due to orbital motion of CTCV 1300 
through the 210 second VIS-arm exposures and rotation- 
ally broadened by a range of velocities (50-200 km s" 1 ). In 
principle, the orbital smearing is a function of orbital phase, 
and thus varies throughout the orbital cycle. We use a single 
value of 46 km s~\ which is the average value of the smear- 
ing across an orbital cycle. We find that changing this to the 
maximum and minimum possible values of orbital smearing 
required, that is the smearing at conjunction and quadra- 
ture, alters the final value of vsini obtained by 3 km s _1 . 
This uncertainty is added in quadrature to the uncertainty 
calculated below. 

The value of v sin i was obtained via an optimal sub- 
traction routine, which subtracts a constant times the nor- 
malised, broadened template spectrum from the normalised, 
orbitally corrected CV spectrum. This constant is adjusted 
to minimise the residual scatter between the spectra. The 
scatter is measured by carrying out the subtraction and 
then computing x 2 between the residual spectrum and a 
smoothed version of itself. By finding the value of ro- 
tational broadening that minimises x 2 , we can obtain a 
value of v sin i and the s pectral type of the secondary sta r 
jDhillon fc Marsbl 1 19931 ; iMarsh. Robinson fc Wood] Il994l ). 
This value of v sin i should then be corrected for the intrin- 
sic rotational velocity of the template star. Unfortunately 
a wide range of spectral-types were not available, and so 
we are unable to deduce the spectral-type of the secondary 
using this technique. 

The value of v sin i obtained using this method was 
found to vary depending on the spectral-type template used 
and the wavelength region selected for optimal subtraction. 
We attempted to include as much of the continuum as pos- 
sible around the Na I doublet, while trying to avoid tel- 
luric regions. We used a wavelength range of 8080-8106, 
8125-8206, 8226-8245 and 8264-8285 A, a limb-darkening 
coefficient of 0.5 and smoothing Gaussian of FWHM = 
15 km s _1 , which were found to give the lowest values of 
X 2 - Th e limb-darkening coefficien t is highly uncertain, al- 
though [CopEerwheat^etjil] l|201ll ) have shown that altering 
the limb-darkening coefficient has little effect on the value 
of v sin i obtained. We plot the values of x 2 versus v sin i for 
both spectral-type templates in Fig. [8] Using the M2 tem- 
plate, we obtain a value of vsini = 129 ± 3 km s _1 , while 
the M5 template yields a value of v sini = 125 ± 4 km s -1 . 
The uncertainties on these values come from the formal error 
estimation of Ax 2 ~ ±1. This does not attempt to include 
systematic errors. Because of the lack of available templates, 
we estimated the spectral type of the secondary sta r using 
the empirical donor sequence of lKnigge et al.l (|201ll ). For a 
system with an orbital period of 128.07 minutes, we expect 
a secondary with spectral type of M4.3. We adopt a spectral 
type of M4.5 ± 0.5. We interpolate between the two values 
of v sin i above to arrive at a final value of v sin i = 125 ± 7 
km s" 1 . This error takes into account both the uncertainty 
on finding our minimum v sin i for each template (± 3-4 km 
s _1 for each template), the uncertainty from averaging the 
orbital smearing (± 3 km s _1 ), and the uncertainty in spec- 
tral type (±0.5 spectral types). 
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Figure 8. % vs usini from the optimal subtraction technique. 



4.6 System Parameters 

Using the values of K2 = 378 ± 6 km s -1 and vsmi = 125 



± 7 km s~ found in Sections 14.41 and 14.51 in conjunction 
with the orbital period and a measurement of the eclipse 
full width at half depth (A^i/2), we can calculate accurate 
system parameters for CTCV 1300. 

The be st measurement of t he orbital period, P OT b, 
comes from ISavourv et al. J201ll). who d etermine P or b = 
0.088940717(1) days ISavourv et all l|201ll ) also present six 
light curves of CTCV 1300, from which we determine A<^i/ 2 
= 0.0791(5). 

We use a Mont e Carlo approach similar to Home et al.l 
I 19931). [Smith et all (1 19981) IThoroughgood et all (j200ll ) and 
iThoroughgood et alj ( 20041 ) to calculate the system param- 
eters and their errors. For a given set of K%, v sin i, P or b and 
A(/>i/2, the remaining parameters are calculated as follows. 

R2/C1 can be estimated because the secondary star fills 
its Roche Lobe. R2 is the secondary radius, and the a is 
the binary sepa ration, and so we use Eggleton's formula 
l|Eggleton| [l983). which gives the volume equivalent radius 
of the Roche Lobe to an accuracy of ~1 per cent, which is 
close to the equatorial radius of the secondary star as seen 
during eclipse, 



R2 
a 



0A9q 



2/3 



0.6g 2 / 3 + ln(l + q 1 / 3 



(2) 



The secondary star rotates synchronously with the orbital 
motion, so we can combine K2 and vsmi, to get 



R 2 , . wsini 

— {1 + q ) = ——. 3 

a K2 

This gives us two simultaneous equations that can be solved 
for q and R2/C1. The orbital inclination, i, is fixed by q and 
using geometrical arguements (e.g. lBailevlll979l ). We 
determine the inclination via a binary chop search using an 
accurate model of the Roche Lobe. 

Using Kepler's Third Law, we obtain 



K2 Por. 

2ttG 



Mi sin 3 



(4) 



which using the previously calculated values of q and i yields 
the mass of the primary star, Mi. The mass of the secondary 
star, M2 and radial velocity of the primary, K\, is given by 

_ M _ Ki 

q ~ Mi ~ K 2 

Finally, we can calculate the radius of the secondary star 
using 



(5) 



v sin i 2tt sin i 



(6) 



R2 Porb 

and the binary separation, a, using equations and |S] 

Our Monte Carlo simulation takes 250,000 values of K2, 
A</>i/2, wsini and P or b, treating each as being normally dis- 
tributed about their measured values with standard devi- 
ations equal to the errors on the measurements. We then 
calculate the mass of each component, the inclination of 
the system and the radius of the secondary star as outlined 
above, omitting (K2, «sini, A(j>i/2) triplets that are incon- 
sistent with sini ^1. Each accepted Mi, M2 pair is plotted 
in Fig. [5] and the masses and their errors are computed from 
the mean and standard deviation of the distribution of these 
pairs. We find that Mi = 0.79 ± 0.05 M and Mi = 0.198 ± 
0.029 Mq. Th ese values are f ound to be in good agreement 
with those of I Savoury et al.l (|201ll ) . The values of all sys- 
tem parameters found from the Mont e Carlo simulat i on are 
listed in Table Q] along with those of I Savoury" et all (|201lD 
for direct comparison. 



4.7 Distance 

By finding the apparent magnitude of the secondary star 
from its contribution to the total light during eclipse, and 
by estimating the absolute magnitude, we can calculate the 
distance (d), using the equation; 



51og(d/10) 



mi 



Mi - dAi/W00, 



(7) 



where Ai is the interstellar extinction in magnitudes per 
kpc. We assume the extinction is zero, as t his allows a di- 
rect c omparison to the distance obtained by I Savoury et all 
(2011), who used model generated white dwarf fluxes to es- 
timate the distance without correction for extinction. At 
mid-eclipse (<j> — 0), the apparent magnitude of the sys- 
tem is 17.32 ± 0.02 around the Na I doublet, which is ap- 
proximately the /-band. This value is not corrected for slit 
losses. The secondary star is found to contribute 58 ± 6 per 
cent, which gives an apparent magnitude of mi — 17.91 ± 
0.09. We estimate the absolute magnitude using the empiri- 
cal donor sequence of lKnigge et~aLl (|201lt ). who assume the 
donor is on the main sequence and then correct for bloating 
effects. From this, we take Mi — 10.32 ± 0.14, and obtain a 
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Figure 9. Monte Carlo determination of system parameters for CTCV 1300. Each dot represents an (Mi, M%) pair. Dot-dashed lines 
are lines of constant inclination, the solid curves satisfy the constraints from the radial velocity of the secondary star, K2, and the dashed 
lines satisfy the constraints of the rotational velocity of the secondary star, i;sini. 



Table 1. System parameters for CTCV 1300. 



Parameter 


Measured Values 


Monte Carlo Values 


Savoury et al. (2011) 


Porb (s) 






0.088940717(1) 




0.0791 ± 0.0005 




not stated 


K 2 (km s" 1 ) 


378 ±6 




372.2 ± 2.5 


vs'mi (km s — 1 ) 


125 ±7 




122 ± 10* 


<? 




0.252 ± 0.025 


0.240 ±0.021 


i° 




85.7 ± 1.5 


86.3 ± 1.1 


Mi/M 




0.79 ±0.05 


0.736 ±0.014 


M 2 /M 




0.198 ±0.029 


0.177 ±0.021 


R2/RQ 




0.223 ±0.011 


0.215 ±0.008 


0,/Rq 




0.834 ± 0.020 


0.813 ±0.011 


A'i (km s- 1 ) 




95 ±9 


90 ±8 


Distance (pc) 


330 ± 40 




375 ± 13 



Derived using the values published in table 3 of lSavourv et al ] l|201ll l. 



distance of d = 330 ± 40 pc. T his distance is found to be 
in good agreement with that of ISavourv et all (|201lf ) . who 
obtained a value of d — 375 ± 13 pc. 



5 DISCUSSION 

The system parameters listed i n Table [T] are found to be 



Savoury et all (2011). To- 



201 ll ). this gives us conn- 



in good agr eement with those o f j 
gether with ICopperwheat et al.l 
dence that photometr ic m ass determinations suc h as those of 
iLittlefair et aH (|2008l ) and lSavourv et alj |201l]) are reliable 
across a range of orbital periods, and that the e-q relations 



of lPattersonl l|2005h . lKnigge] l|2006|) and lKnigge et all (|201lD 
are well founded. 

The uncertainties in the system parameters for CTCV 
1300 determined in this paper, and in ISavourv et ail i|201ll ). 
are quite lar ge in comparison to m any of the other systems 
published in ISavourv et al. The reason for the large 

uncertainties in CTCV 1300 in lSavourv et alj (|201ll ) is be- 
cause the eclipse light curves used for model fitting suffer 
from heavy flickering, which causes difficulties in obtaining 
an accurate value for the mass ratio, q. The large uncertain- 
ties in this paper arise because of the interpolation technique 
used to arrive at a value for usini. The error on vsini (± 
7 km s" 1 ) is the dominant source of uncertainty in our final 
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system parameters. In principle, a wider selection of spectral 
type templates would enable us to further constrain vsini, 
and derive the spectral type. 



6 CONCLUSIONS 

We have used time-resolved spectroscopy to determine the 
system parameters for the short period dwarf nova CTCV 
1300. The double-peaked nature of the Balmer and He I lines 
confirms the presence of an accretion disc, while careful anal- 
ysis of the Na I doublet absorption lines at 8189 A reveals 
the radial velocity of the secondary star to be K2 = 378 ± 
6 km s _1 and the rotational velocity of the secondary star 
to be vsini = 125 ± 7 km s . Using these measurements, 
we find Mi = 0.79 ± 0.05 Mq for the white dwarf primary 
and Mi = 0.198 ± 0.029 M© for the M-type secondary star. 
The radius of the secondary star is found to be R2 = 0.223 
± 0.011 R Q . 

The system parameters determined through spectro- 
scopic analysis are found to be in good agreement with those 
previously calculated using photometric techniques. This is 
significant, as our results support the validity and accuracy 
of the purely photometric mass determination technique in 
short period cataclysmic variables. 
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